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§ i ABSTRACT 
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o. 

Using the spatial classification method and the structural parameters esti- 
mated for the Gould Belt (GB) and the local Galactic disk (LGD) from a pre- 
vious paper, we have evaluated spatial membership probabilities for a sample of 
and B stars from the Hipparcos catalogue (ESA 1997) with available proper 
motions and radial velocity data. Thus being able to study the space velocity 
fields of both systems, we conclude that the GB and the LGD present different 
statistical distributions, both in the velocity space as well as in the phase space. 
^ ■ In the light of the existence of both systems we analyze different kinematic as- 

pects in the vicinity of the Sun, and we find the GB to be responsible for the 
highly negative vertex deviation found for the OB stars in the solar neighbor- 
hood. We also find that the GB sensibly alters the estimation of Oort's constants 
in the Galactic disk; thus, in order to establish comparisons with other kinematic 
studies based in older stellar populations, a careful pruning of the GB members 
must be performed. Further analysis of the GB velocity field and the moving 
groups that contribute to it suggest the possibility that the GB can be the result 
of a chance superposition of moving groups. We propose for future investigations 
the dynamical analysis of these moving groups in order to probe the origin of the 
GB. 



Subject headings: (Galaxy:) solar neighborhood — stars: early-type — stars: 
kinematics 
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1. Introduction 

Since the discovery of the GB as a system of bright stars forming an angle of about 
20° with the Galactic disk (Gould 1879; Herschel 1847), many efforts have been devoted to 
unravel its complex structure of stars and the associated interstellar medium. When studies 
of its kinematics began to appear on the second half of the past century, the peculiarities of 
its behavior made the global picture of the GB even more puzzling. 

The most striking discovery was the expansion motion of the stellar component (Bon- 
neau 1964; Lesh 1968; Fricke & Tsioumis 1975; Frogel & Stothers 1977; Westin 1985; Comeron 
et al. 1994), which cannot be satisfactorily explained by a single explosive event (Lesh 1968; 
Moreno et al. 1999), making it difficult to trace back the origin of the GB by reversing this 
movement. Models considering Galactic density waves generated by perturbations from the 
spiral arms cannot account for this expansion either (Westin 1985; Comeron & Torra 1991). 
Also, the shell of gas associated with the GB present evidences of an expansion movement, 
according to the works of Lindblad (1967), Olano (1982), Elmegreen (1982) and Moreno et 
al. (1999). 

Because the age of the GB (between 20 and 90 Myr, see Torra et al. (2000) for a 
detailed discussion) is a considerable fraction of the period of vertical oscillation of the stars 
over the Galactic plane under the influence of the gravitational potential of the Galaxy, the 
spatial coherence of the GB must be corresponded by a kinematical coherence that prevents 
the dissolution of the structure among the Galactic plane. Comeron (1999) and Perrot & 
Grenier (2003) find a global oscillation movement of the GB around an axis different from 
the spatial line of nodes where the GB and the Galactic plane cross each other. 

All this leads Elmegreen et al. (2000) to include the GB as part of the hierarchy of stellar 
complexes within the Milky Way, as a second level structure subordinated to the local arm 
(the Orion-Cygnus spur). Thus, the GB would be the star formation complex closest to us, 
so the study of its global properties will contribute to shed light on the possible origin and 
evolution of these complexes. 

In a previous article (Elias et al. 2006, from now on referred to as Paper I) we had 
developed a three-dimensional spatial classification method to separate the GB stars from 
the LGD stars. Having obtained the structural parameters of both systems exclusively by 
spatial considerations, we now shall classify a sample of Hipparcos (ESA 1997) OB stars 
with space velocities, in order to compare the kinematics of both structures. 

Thus, in Section 2 we first build a sample of 0-B6 Hipparcos stars with proper motions 
and radial velocity data, that we then analyze in order to study the kinematic properties of 
the young Galactic disk. We begin (Section 3) by the identification of the moving groups 
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in the velocity field, and then (Section 3.1) we separate the GB from the LGD stars, which 
allows us to decide about the membership to either system of the detected moving groups. 
The next step to enhance our analysis of the young Disk velocity field is the elimination 
of the systematic effects on the velocities by the solar motion and the Galactic differential 
rotation. We obtain the residual velocities, the analysis of which also yields a different 
kinematic behavior of the GB and the LGD (Section 4). The study of the velocity ellipsoids 
confirms such difference (Section 4.1). The effect of the GB on the determination of the Oort 
constants is addressed in Section 5. Finally, some conclusions are exposed in Section 6. 



2. Star sample 

We select a sample of 1156 stars of spectral types 0-B6 and luminosity classes III, 
IV and V from the Hipparcos catalogue (ESA 1997). Photometric data from the Hauck 
& Mermilliod (1998) catalogue, as well as radial velocity data from the Barbier-Brossat & 
Figon (2000) and Grenier et al. (1999) catalogues have been added when available. Thus, 
the compilation includes: 

• HIP, Hipparcos identifier number 

• Spectral type 

• V, Johnson visual magnitude 

• Trigonometric parallax (mas) 

• Standard error in trigonometric parallax (mas) 

• Right ascension for the epoch J 199 1.5 in ICRS (degrees) 

• Declination for the epoch J1991.5 in ICRS (degrees) 

• Proper motion in right ascension, ji a cos5 (mas yr _1 ) 

• Standard error in proper motion in right ascension (mas yr _1 ) 

• Proper motion in declination, [is (mas yr _1 ) 

• Standard error in proper motion in declination (mas yr _1 ) 

• Radial velocity (km s -1 ) from Barbier-Brossat & Figon (2000) 

• Quality of radial velocity from Barbier-Brossat & Figon (2000) 
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• Radial velocity (km s -1 ) from Grenier et al. (1999) 

• Error in radial velocity (km s" 1 ) from Grenier et al. (1999) 

• uvbyf3 photometry data from Hauck & Mermilliod (1998) 

For the distance estimation we have used Hipparcos trigonometric parallaxes only if the 
relative error is lower or equal to a 10%. Otherwise, uvbyf3 Stromgren photometry has been 
used to estimate the distance through the Balona & Shobbrook (1987) My(f3) calibration. 
If no data were available, spectro-photometric distances from the apparent visual magnitude 
V and the Schmidt-Kaler (1982) calibration for the spectral types have been chosen. We 
have compared the three different distance estimations, looking for any possible systematic 
biases among them. With this purpose, we have selected 950 stars from the initial catalog 
with Stromgren photometry data. 

Several authors (Arenou & Luri 1999; Maiz- Apellaniz 2001; Maiz Apellaniz 2005; Schroder 
et al. 2004), studying the distance calibration comparisons with the distances obtained from 
trigonometric parallaxes, have dealt with the problem of analyzing the biases that the sam- 
ple selection effects introduce. The estimation of these biases is very complex because it 
depends, among other variables, on the spatial distribution of the sample. Maiz Apellaniz 
(2005) demonstrates that the real distance probability distribution for individual stars will 
always be ill-behaved when that distance tends to infinite and a constant underlying spatial 
distribution has been assumed for the sample. This idea had already been suggested by 
Arenou & Luri (1999), who propose that -in order to avoid any truncation biases- for com- 
parison purposes it should be used a sample not selected by any limit in the relative error of 
the trigonometric parallax, including also the negative parallaxes. We use this methodology 
for the comparison between trigonometric, photometric and spectroscopic parallaxes, using 
the complete sample with spectral types up to B6. 

The result of the analysis demonstrates that the Stromgren photometric parallaxes (iistr) 
are very similar to the Hipparcos trigonometric parallaxes -as Kaltcheva & Knude (1998) 
had already demonstrated-, but the former present, in comparison with the spectroscopic 
parallaxes (tisk, derived from the Schmidt-Kaler calibration), a functional relationship in the 
form of 7istr/^SK ~ 1.21. Nevertheless, a recent study for the O stars (Maiz Apellaniz 2005) 
demonstrates that the Hipparcos trigonometric parallaxes and the spectroscopic parallaxes 
are very similar for this spectral type. Where does this difference that we find in their values 
come from? We have to consider that our spectroscopic calibration is based on three steps: 
1.- Spectral classification; 2.- Calibration of the absolute magnitude for each spectral type; 
3.- Evaluation of the reddening from the intrinsic color values. 
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The spectral classification of the Hipparcos catalog comes from different sources and 
thus is far from being uniform. On the other hand, the catalog of O stars (Maiz-Apellaniz 
et al. 2004) used by Mafz Apellaniz (2005) for the comparison of the parallaxes was derived 
from a group of homogeneous and precise stellar spectra. Thus, it must not be inferred from 
our result that the Schmidt-Kaler calibration has systematic biases, but that the spectro- 
scopic parallaxes of the 0-B6 stars from the Hipparcos catalog -according to the spectral 
classification within the catalog- show a systematic error when we compare them with the 
Stromgren photometric parallaxes (which, as we said, bear no significative difference with 
the Hipparcos trigonometric parallaxes). 

In this work, we do not intend to perform an exhaustive study of the problems that 
the diverse methodologies of obtaining distance calibrations arise. We have just unified for 
working reason the different distance calibrations used for our sample. For that, we have 
tied the spectroscopic parallaxes to those derived from Stromgren photometry, applying a 
correction of 21% to the former. We want to stress that the comparison between the different 
distance calibrations has been made strictly through their respective parallaxes. 

We have also chosen the Barbier-Brossat & Figon (2000) radial velocity and quality 
data when available; otherwise, radial velocities and errors from Grenier et al. (1999) have 
been used. Finally, a distance limit of 1 kpc has been imposed, thus reducing the sample to 
881 stars. 

While the Hipparcos catalogue is complete down to V ~ 7.9, and to V ~ 7.5 for 0-B6 
stars, in Figure 1 we can see through a histogram that our sample is complete only down 
to a magnitude of V ~ 6.5. This is caused by the necessity of having radial velocity data 
available for the stars in our sample, as it is also observed in Fernandez (2005); in that work, 
the completeness of the sample of O and B stars from the Hipparcos catalogue falls down 
from V ~ 7.9 to V ~ 6 when the stars without radial velocity data are removed. 

3. Identification of moving groups in the sample and their membership 

We have calculated the space velocities from the proper motions and the radial velocities 
(Johnson & Soderblom 1987) for the stars in our sample. Their density field is represented 
in the three contour plots of Figure 2. We appreciate in this image that the velocity field is 
dominated by several maxima which may correspond to associations of stars (not necessarily 
bound) with a small velocity dispersion, known as moving groups (e.g. Proctor 1869; Eggen 
1963). 

The first and most prominent maximum, located around (U,V,W) = (—6.5,-19,-7) 
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km s -1 , is certainly associated with the Pleiades moving group. The exact situation of the 
peak may differ slightly from the estimation given by other authors, but we must consider 
that the moving group always appears as a maximum in the velocity space with a certain 
width. For instance, Chen et al. (1997), working with a sample of B, A and F stars from 
the Hipparcos Input Catalogue (Turon et al. 1992) and uvbyf3 photometry, find that the 
maximum is located at (U, V, W) = (—10, —19, —8.1) km s _1 , the standard deviation for the 
three components being, respectively, (7.9,8.6,5.8) km s _1 . 

The second most prominent peak in Figure 2 is located around (U, V, W) = (—17, —11, —5) 
km s -1 . We have identified it as the moving group related to the supercluster IC 2391. Its po- 
sition in velocity space is estimated by Chen et al. (1997) at (U, V, W) = (—15.9, —13.1, —4.5) 
km s -1 , the standard deviation being (<Ju, ay, a w ) = (4.1,6.2,3.0) km s _1 . 

The third maximum, centered around (U,V) = (—11,-8.5) km s -1 , is more diffuse 
and difficult to identify. We have decided to follow the criteria of Asiain et al. (1999), who 
-based on the studies by Comeron (1992)- rule out the possibility of linking this region in 
the velocity space with the Coma Berenices cluster, in favor of a probable bond with the 
Cassiopeia- Taurus association, of (U,V) = (—9.9,-6.1) km s _1 . This association has an 
age between 50 Myr -as the probable expansion age estimated by Blaauw (1956) from a 
sample of stars of spectral type B5 or earlier- and 90 ± 10 Myr, as determined from the 
Lithium depletion boundary by Stauffer et al. (1999). In the detailed study of the nearby 
OB associations by de Zeeuw et al. (1999) the authors find a physical relation between the 
Cas-Tau group and the a Persei cluster. The main sequence turn-off age for this cluster is 
about 50 Myr (Meynet et al. 1993), so the age of both Cas-Tau and a Persei could be the 
same (e.g. Brown 2002). Note that in Paper I we found that a Persei is located well within 
the spatial boundaries of the GB. 



3.1. Classification of the sample 

In order to observe how these moving groups contribute to the two systems in study, 
we classify the sample and determine the membership of the stars to either the GB or the 
LGD. In Paper I we obtained several estimations of the parameters that characterize both 
systems in our model. We shall work with the solution for the full 0-B6 sample of Paper I 
with an exponential model for the stellar distribution in the direction perpendicular to the 
mean planes. 

Thus, using that estimation as the true value of the GB and the LGD mean planes, we 
simply apply our separation algorithm to our current star sample with kinematic data. No 
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iteration in order to re-evaluate the mean planes is performed, we just assign membership 
probabilities to the stars according to the planes estimated in Paper I. Also, spatial outliers 
are eliminated according to the procedure explained in that paper; the remaining sample has 
776 stars. We obtain a separation between the GB and the LGD based exclusively in the 
spatial position of these stars. Yet we can see in Figure 3 how a difference in their velocity 
fields is obtained as a result. 

The most striking difference is that in the UV projection (top panels of Figure 3) 
the three moving groups that we had found in the full sample (top panel of Figure 2) 
now distinctly belong either to the GB or the LGD. The two maxima associated with the 
Pleiades and IC 2391 appear only in the GB velocity field (Figure 3, top left panel), while 
the Cassiopeia- Taurus peak remains only visible in the LGD field (Figure 3, top right panel). 
This is not surprising, if we consider that the Pleiades moving group is spatially related to 
the Sco-Cen association, which is one of the main components of the GB (Moreno et al. 
1999). Also, we know that IC 2391 is a young cluster, its age being about 30 Myr (Stauffer 
et al. 1997). In a recent study of the tangential velocities, Piskunov et al. (2006) conclude 
that its kinematic probability of belonging to the GB is a 73%. Note that we have arrived to 
a similar conclusion by a process based solely on the spatial position of the stars, and thus 
independent of the result in the cited paper. 

We must also note that a late-type population of young stars has been associated 
to both the Pleiades and IC 2391 by Montes et al. (2001). In that paper, these moving 
groups are described as centered around (U,V,W) = (—11.6,-21,-11.4) and (U,V,W) = 
(—20.6,-15.7,-9.1) km s -1 , respectively, with a dispersion of about 8 km s _1 around the 
central positions. Not surprisingly, a late-type component of stars of about 30-80 Myr of age 
had already been associated to the GB disk structure by Guillout et al. (1998) studying the 
X-ray sources in the ROSAT All-Sky Survey. 

Finally, we also observe that in the UW projection (Figure 3, middle right panel) it 
clearly rises a new maximum that has a correspondence with a small protuberance around 
U = 11 km s -1 in the UV projection (Figure 3, top right panel). It was also present, although 
weak, in the velocity density field of the full sample (Figure 2, top and middle panels). We 
haven't found an exact correspondence to this possible moving group among the structures 
in the solar neighborhood, but the positive value of the U component makes us think that it 
may be related to the Sirius supercluster (eg. Eggen 1996, Dehnen 1998, Asiain et al. 1999). 
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4. Analysis of the residual velocities 

If we want to improve our study of the kinematic differences between the GB and the 
LGD, we must work with residual velocities. This way, the systematic effects of the Galactic 
kinematics will not interfere with our comparison. Thus, we now correct the space velocities 
of the stars in our sample from solar motion (using the classical estimation of Delhaye (1965), 
([/©, V , W Q ) = (9, 12, 7) km s" 1 ) and from Galactic differential rotation -using the values 
given by Kerr & Lynden-Bell (1986). 

In order to refine the analysis of the velocity distributions, we must also eliminate 
outliers of kinematic nature that may be present in the sample. That is, we eliminate the 
stars located in the regions of a very low density in the velocity space for our sample under 
study. We achieve this by running the OUTKER algorithm (Cabrera-Caho & Alfaro 1985), 
which reduces the sample to a final number of 752 stars. We then separate again the sample 
into GB and LGD members. Now we want to compare both distributions in a N-dimensional 
space, where N=6 (phase space) or N=3 (velocity space). In order to achieve this, we employ 
a multidimensional, non-parametric two-samples test: the Cramer test (Baringhaus & Franz 
2004). The test statistic is the difference of the sum of all the Euclidean interpoint distances 
between the random variables from the two different samples: 

i,j i,j=l i,j=l / 

(1) 

where Xi and Y i are the point vectors of each sample member, and m,n are the 
respective size of the samples. is a kernel function; for this particular case we have 
used the Cramer kernel implemented by Franz (2004) in the R statistical environment (R 
Development Core Team 2005). 

Thus, an analysis of the three-dimensional velocity space, (U, V, W), yields that the GB 
and the LGD distributions are different with a 99% of confidence. Similarly, the Cramer test 
for the six dimensions of the phase space, (X, Y, Z, U, V, W), rejects the possibility that the 
GB and the LGD distributions are the same with a confidence of a 99%. 

Thus we confirm that the GB and the LGD are two different stellar systems in the sense 
that they show a clear statistical separation between their distributions in the phase space. 
It is unavoidable to take into consideration the contribution of the GB when the young Disk 
is under study, because, as we have demonstrated, the velocities of the OB stars in the solar 
neighborhood are not statistically independent of their spatial positions. 
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4.1. Velocity ellipsoids 

A more intuitive visualization of the differences between the GB and the LGD residual 
velocity distributions are the contour density plots of the three different velocity planes (see 
Figure 4). Their disparate shapes and orientations already show that the velocity ellipsoids 
clearly reflect the distinct kinematic behavior of both stellar systems. 

We have estimated the main geometrical parameters of the velocity ellipsoids for the 
whole sample and for the GB and the LGD members separately. The results are displayed 
on table 1. Two main results arise from this analysis: 

f . The vertex deviation for the GB is highly negative, while that for the LGD is positive. 

2. There is the suggestion that the third axis of the velocity ellipsoid for the GB is tilted 
in respect to the Galactic plane. 

The estimation of the vertex deviation in the solar neighborhood from different star 
samples has produced different values depending on the nature of the sample and on the 
kinematic variables used in the calculation (see Moreno et al. (1999) for a compilation of 
previous results). In brief, the general conclusion has been that the vertex deviation becomes 
more negative as the star sample is younger. In fact, one of the classic estimates based 
on space velocities for OB stars (Filin 1957) yielded a value close to l v = —50° for the 
Galactic disk. For years, this result has remained a puzzling issue that has been given 
several and varied explanations. Most of them can be classified into two classical types: 
nature or nurture, we could say. Some authors claim that these young stars that formed 
from a molecular cloud show the same kinematics as the parent cloud at the time of the star 
formation. This way, the initial velocity and later expansion define the velocity ellipsoid of 
the present stellar system. Other authors, though, defend that the effect of different singular 
and punctual events (such as the passing through a spiral arm) could also be the cause of 
the peculiar velocity ellipsoid observed in the young stellar component. 

What we conclude from our analysis is that the classic problem of the negative vertex 
deviation for young stars in the solar neighborhood is a consequence of the presence of the 
GB. If we eliminate the stars that belong to the GB, the remaining sample of only LGD 
stars present a positive vertex deviation (l^ GD = 18°). 

Moreno et al. (1999), working with a sample of dwarf 0-B5.5 stars members of the GB, 
found that the negative vertex deviation (l v ~ —64°) was caused by the Pleiades moving 
group. Once removed this group, they obtained a positive vertex deviation (l v = 22°) for 
the remaining stars in the GB. We now demonstrate the the OB stars of the LGD also 
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have a positive vertex deviation, similar to that found by Moreno et al. (1999) for the GB 
stars after eliminating the ones belonging to the Pleiades moving group. So when we work 
with samples (either of GB or LGD stars) in which a single moving group clearly dominates 
over the others the vertex deviation is positive and close to l v ~ 20°, which is the value 
expected from the dynamic equations of the system for this age group. The negative vertex 
deviation of the GB seems to be originated by the relative position of the centroids of two 
moving groups rather than by the distribution of the residual velocities as a whole. We must 
note that the analysis of the young stars -specially those belonging to the GB- based on the 
Schwarzschild distribution from which the velocity ellipsoid is tailored, is not the best suited 
to describe the reality of the system. The velocity field is dominated, as we have seen, by 
the presence of moving groups, and thus it is far from the hypothesis of a homogeneous and 
stationary system. The problem thus lies in that the classic analysis of the velocity ellipsoid 
is applied to a set of moving groups. Already Mihalas & Binney (1981) pointed out that 
the cause of the vertex deviation was the existence of moving groups in the Galactic velocity 
field. But even though the velocity ellipsoid does not strictly correspond to a physic reality 
in our case, that doesn't invalidate our result, that is, the fundamental contribution of the 
GB to the negative vertex deviation for the O and B stars. 

Thus, the different values of the vertex deviation found in the literature can be explained 
according to the different proportions of GB stars present in the respective samples. This 
translates the question about the origin of the negative vertex deviation to the investigation 
of the origin of the moving groups. Although the latter is out of the scope of this paper, 
we want to stress that the more we study in detail the nature of the GB, the more we find 
indications that we must probe both its nature and origin as a set of moving groups. 

Another striking result is the inclination of the GB's velocity ellipsoid of about 10° ± 
12° with respect to the UV plane, although this value does not have a great statistical 
significance. While the third axis (W) of the LGD ellipsoid merely shows an inclination of 
2°, the GB's ellipsoid inclination (10°) reminds that of the spatial system with respect to 
the Galactic plane (icB ~ 14° as we found in Paper I). Some authors have pointed out that 
the GB could be oscillating as a whole around the Galactic plane (Comeron 1999; Perrot 
& Grenier 2003). But our result does not discard a different interpretation in terms of a 
juxtaposition of moving groups: this links back to the idea that a descriptive parameter of 
the GB -such as this inclination of the minor axis of the ellipsoid- can be interpreted in terms 
of the relative position (in the velocity space) of two moving groups. 
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5. Estimation of the Oort constants 

Although our sample hasn't been ideally compiled with the intention of calculating the 
Oort constants, it is worth performing a basic estimation to obtain information about the 
GB's contribution to their value. Thus, within the axisymmetric approximation of Oort's 
model in first order, and including a K expansion term, we solve the condition equations 
(Smart 1968; Clube 1972; Frogel & Stothers 1977) for the radial and the tangential velocities: 

v r = u cos I cos b + v o sin I cos b + w sin b + Ar sin 21 cos 2 b + K (2) 

v t = —-u sin I + v cos I + Ar cos 21 cos b + Br cos b (3a) 
Vb = —Mo cos I sin b — vq sin I sin b + wo cos b — Ar sin 21 cos b sin b (3b) 

where i> r is the radial velocity; vi = 4.74057 and v b = 4.74057 « 6 r are the tangential 
velocities in Galactic longitude (/) and latitude (6), being and u& the respective proper 
motions and r the heliocentric distance; and («o, t>o, m>o) = — (U q ,V q ,Wq), the reflex of the 
solar motion. 

In principle, it would only make sense to estimate Oort's constants for the LGD, elim- 
inating the GB members from the sample. As we have seen, from the study of the moving 
groups and the velocity ellipsoid, the velocity distribution of the GB clearly deviates from 
the axisymmetric hypothesis; hence, the Oort constants wouldn't correspond to a physical 
reality in this case. Yet we have solved the equations in order to establish a comparison 
and thus observe the effects that the presence of the GB introduces in the kinematics of the 
LGD. 

The solutions for the proper motions and for the radial velocity alone are listed on 
table 2. The solar motion for the LGD proper motions solution is in good agreement with 
the IAU standard, (U e , V e , W & ) = (9, 12, 7) km s~ l (Kerr & Lynden-Bell 1986). But the 
Oort constants differ noticeably from the IAU recommended values of A = 14.4 ±1.2 km 
s _1 kpc -1 and B = -12.0 ± 2.8 km s _1 kpc -1 , and instead we have a flat rotation curve with 
A = —B = 16 km s _1 kpc -1 . Yet the value provided by Kerr & Lynden-Bell (1986) is a 
mean over the results obtained by several authors, so it is more revealing to compare our 
estimation with that of a study similar to ours. 

Frogel & Stothers (1977), working with a sample of 0-B5 stars, separate the GB from 
the LGD, and considering a fixed value of A = 15 km s -1 kpc -1 , they estimate the value 
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Table 1. Velocity ellipsoids 



<7 


I 


6 


Axis (km s _1 ) 


(°) 


(°) 



Full sample: 

£/' 10.7 ±0.4 

V 10.2 ±0.4 
W 6.9 ±0.6 

Gould Belt: 

U' 9.8 ±0.3 

V 8.9 ±0.3 
W 7.1 ±0.8 

Local Galactic disk: 
U' 12.2 ±0.6 

V 10.2 ±0.5 
W 7± 1 



-5 ±29 6 ±4 
85 ±29 -2 ±4 
153 ± 115 84 ±3 



-47 ±22 -10 ±8 
43 ± 22 1 ± 12 
141 ±92 80 ± 12 



18 ±11 -2 ±4 
-72 ±11 0±5 
27 ± 66 88 ± 5 



Note. - U', V and W' represent the 
principal axes of the ellipsoids, which devi- 
ate from the (U, V, W) reference frame. The 
errors were estimated by bootstrap. 



Table 2. Oort constants 



Uq Vq Wq a b k 

Sample (km s _1 ) (km s _1 ) (km s _1 ) (km s _1 kpc -1 ) (km s _1 kpc -1 ) (km s _1 ) 



Solution from proper motions: 

FS 9.8 ±0.4 13.0 ±0.6 6.6 ± 0.3 

GB 9.9 ±0.5 13.0 ±0.6 6.7 ± 0.3 

LGD 9.4 ±0.8 12.6 ±0.9 6.3 ± 0.3 

Solution from radial velocity: 
FS 8.7 ±0.8 15.0 ±0.6 11 ±2 

GB 10 ±2 15 ±1 8 ±4 

LGD 9±1 14 ±1 8±5 



14 ±1 -18 ±1 

11 ±2 -20 ±1 

16 ±2 -16 ±1 

13 ±2 - -0.6 ±0.4 

9±2 - 0.4 ±0.6 

16 ±3 - -2±1 



Note. — FS is the Full Sample 
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of B for both systems, as well as for the whole unclassified sample. For the latter they find 
a flat rotation curve of A = —B = 15 km s _1 kpc -1 ; for the LGD, B = -12 km s _1 kpc -1 ; 
and for the GB, B = -19 km s -1 kpc -1 , a value similar to our estimation. But regardless 
of the numerical results, we can conclude too that the global kinematics of the GB differs 
greatly from the kinematics of the LGD, and that the presence of this stellar system affects 
the estimation of the parameters describing the kinematics of the solar neighborhood. 

On the other hand, the K term is generally inferior in absolute value to the estimation 
by Frogel & Stothers (1977). We have checked that when solving the equations without 
this term, the results for the Oort constants do not change significatively. Thus we do not 
consider it wise to conclude anything about the possible expansion movements of the system 
from these results. 

The detailed study of the local irregularities in the kinematics of the young stars by 
Torra et al. (2000), from a sample of O and B Hipparcos stars, reveals that for heliocentric 
distances lower than 600 pc and age groups under 90 Myr (i.e., for a sample with a high 
proportion of GB members), the value of the Oort constant B is much greater in absolute 
value than expected for the LGD: -13.6 ± 2.0 < B < -20.7 ± 1.4 km s -1 kpc -1 . On the 
other hand, they find a quite small value of the Oort constant A: 10.5 ±2.1 < A < 5.7 ± 
1.4 km s -1 kpc -1 . Both estimations are perfectly compatible with our result of (A, B) = (11 
± 2, -20 ± 2) km s -1 kpc -1 for the GB, which is undoubtedly contaminating their sample. 

Some other recent studies find values of A and B in good agreement with ours. Uemura 
et al. (2000), for a sample of 0-B5 Hipparcos stars, estimate that A = 14.0 ± 0.7 km 
s -1 kpc -1 and B = -15.8 ± 0.7 km s -1 kpc -1 ; Zhu (2000), working with 0-B5 stars with 
Hipparcos proper motions, finds that A — 16 ± 1 km s -1 kpc -1 and B = -15.6 ± 0.8 km s -1 
kpc -1 ; Oiling & Dehnen (2003), for a sample of young main sequence stars find that (A, B) 
= (9.6, -11.6) km s -1 kpc -1 , while for a sample of red giants, with no significant contribution 
from the GB, the Disk rotation curve is almost flat, (A, B) = (15.9, -16.9) km s -1 kpc -1 , in 
very good agreement with our results for the OB stars of the LGD. 

Our estimation of the circular rotation speed, considering that the Sun's distance to 
the Galactic center is 8.5 kpc (Kerr & Lynden-Bell 1986), is O = 272 km s -1 for the LGD 
from proper motions. It is a high value, but in agreement with recent studies of 0-B5 star 
samples from the Hipparcos catalog by Miyamoto & Zhu (1998), O = 268.7 ± 11.9 km 
s -1 ; Uemura et al. (2000), O = 255.52 ± 8.33 km s -1 ; or Branham (2002), O = 258.7 ± 
34.29 km s _1 . A similar result (O = 270 km s -1 ) is also reached by Mendez et al. (2000) 
using data from the Southern Proper Motion catalog (Platais et al. 1998); thus, their result 
is obtained from measures of the proper motions independent from those of Hipparcos. 
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6. Conclusions 

In Paper I we concluded that in the young Galactic disk there are two spatially different 
systems, the GB and the LGD. Now we have classified only by spatial criteria a sample 
of 0-B6 stars from the Hipparcos catalogue, and found that the GB's global kinematics is 
essentially different from the kinematics of the LGD. Not only that, multidimensional two- 
samples tests prove that the GB and the LGD are separated systems in the phase space. This 
does not necessarily imply that the GB is a coherent structure born from a single source such 
as a great molecular cloud. We just note that in its present state, the GB is a local system 
(whose size is well represented by our sample) showing very different kinematic properties 
than the larger system in which it is embedded, the Galactic disk. The fact that the GB is 
mainly composed of certain moving groups challenges the idea of this system coming from a 
single origin, and rises the question of whether we are witnessing a coherent structure with 
physical entity or we are just observing a transitory picture of several smaller systems with 
no common origin at all. The answer to this question can only be sought in the dynamical 
study of the moving groups that form the GB (which would require a deep knowledge of the 
Galactic potential and its asymmetries), in order to trace back their trajectories in the past 
and discover whether they come from a single protostellar cloud or they are the transitory 
result of some dynamical traps. 

We have also proved that the classic problem of the negative vertex deviation of young 
stars in the solar neighborhood is caused by the contribution of the GB. This effect disappears 
once this stellar system is removed from the sample, leaving the LGD -defined by its spatial 
distribution- as the only remaining structure. 

Finally, we have observed how the presence of what is called the GB introduces distur- 
bances in the estimation of the Oort constants that describe the kinematics of the young 
Galactic disk, making it necessary to discard its contribution by identifying and removing 
the GB members. Once the young Disk is pruned by eliminating the GB members, we esti- 
mate a flat rotation curve with a local velocity very close to the values calculated by several 
authors in the last decade for a wide range of ages. 

Thus, although we find that the kinematic analysis is not enough to decipher the origin 
of the GB, it is indeed fundamental to characterize the complexity of the young Galactic disk 
and to better understand the different moving groups that form the bulk of the GB stellar 
component. A comprehensive explanation of the origin of the GB will require a dynamical 
analysis of these moving groups. 
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Fig. 1. — Histograms of the V magnitude in the Johnson system for the 0-B6 stars in the 
Hipparcos catalogue (ESA 1997) and for our star sample. 
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Fig. 2. — Distribution of space velocities for the star sample. 
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Fig. 3. — Distribution of space velocities for the GB (left) and the LGD (right) stars. 
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Fig. 4. — Contour density plots of the residual velocities for the GB (left) and the LGD 
(right) stars. 
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